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^ : Abstract 
cn 
o 

' The low-mass dielectron signal from heavy-ion collisions at the CERN-SPS 

0^ ' reported by the CERES collaboration is in excess of estimated hadronic de- 



cays suggestive of possible contribution from two-pion annihilation or other 
hadronic reactions. In the absence of dramatic medium modifications, an- 
nihilation alone is unable to account for the data. We explore the role of 



^ ' pion plus resonance scattering [irp ai(1260) — > ne'^e ] which has favorable 



kinematics to populate masses between 2m.j^ and rup. While it seems to ac- 
count for some of the remaining excess beyond annihilation, it fails to allow 
quantitative interpretation of data. 
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Dynamics of high-energy heavy- ion colhsions can be probed, at least in principle, by 
photons and dileptons since they are produced continuously from the first moments of nuclear 
contact throughout the entire evolution and interact only weakly through electromagnetic 
coupling to the hadronic medium. Their mean free paths are many times greater than 
typical reaction-zone sizes allowing them to bring out valuable information on the earlier 
stages. Since their utility for heavy- ion collisions was first proposed many theoretical 
calculations have been done predicting orders of magnitude for production, exploring 
ideas about phase transition to quark matter and possible medium modifications to hadronic 
phenomena. Recently an intriguing report from CERN brought forth results for dielectron 
signals observed in heavy-ion induced reactions as compared with proton-induced results. 
While the p-l-Be and p-l-Au data at 450 GeV were not inconsistent with estimations from 
hadronic decays, the S-I-Au results showed an enhancement by a factor of 5.0± 2.7 over 
hadronic decay contributions when integrated over mass from 0.2 to 1.5 GeV. 

The natural conclusion drawn from these theoretical estimates and experimental data was 
that the excess is most probably due to tttt annihilation in a heated hadronic environment. 
But since the excess is pronounced for masses near 500 MeV and lower, this is not obviously 
annihilation. If not, then what? Perhaps we are seeing the rho mass shifting downward by 
several hundred MeV due to the presence of the medium ; or perhaps we are seeing effects 
of modified pion dispersion relations |@]. Both phenomena would be spectacular if only we 
could unequivocally establish their influences which at this stage are matters of discussion. 
In fact, prediction for direction and magnitude of mass shifts at finite energy density are not 
yet settled as the literature bears opposing results P,p!0[|. Taking the point of view to use 



free space masses and widths (and including a collision broadening for p), we show in this 
letter that the "excess" dielectron signal is due in part due to tttt annihilation, but could 
have important contributions from pion-resonance irp ne^e~ scattering. 

Consider a central collision of sulfur on gold at 200 GeV/u which will produce a "fire 
cylinder" of transverse radius Rt equal to that of the sulfur nucleus. The question of initial 
energy densities and QCD phase transition are quite interesting but will not be discussed 



in any detail here. Only to note that assuming isentropic expansion and relating observed 
rapidity densities to initial times and temperatures 0, Tj is quite high — of the order 350 
MeV. Instead, if either the mixed phase (of quark and hadron coexistence) or pure plasma 
phase is reached, the additional degrees of freedom will effectively keep the temperature 
lower. This is the picture envisioned so that Tj is of the order 180 MeV. Simulations of 
relativistic transport support this picture since they are able to reproduce transverse 
mass spectra for both pions and protons with an initial temperature of ~ 185 MeV and a 
finite chemical potential = 135 MeV. Being only somewhat above the critical value of 
temperature obtained from lattice calculations fill for entering the plasma phase, we neglect 
contributions from quark processes. The system expands and cools very quickly until the 
mixed phase is reached when hadronic contributions begin thermally producing dielectrons 
at Tc for a few fm/c until the matter is found to be in a complete hadronic phase. It cools 
and expands further until collisions can no longer support the nearly isotropic momentum 
distributions and freezes out. 

The observed spectrum of dielectrons will consist of several components, each coming 
from different mechanisms and possibly different stages of the evolution. Below 150 MeV 
mass a nonthermal component of vr^ and rj Dalitz dominates and the observed yields are 
consistent with estimates |^ . Hadronic radiative decays will contribute to higher masses and 
the largest seems to be ^ 7r°e+e~ []6|,p!2|, although there are recent suggestions that partial 



U{1)a restoration could modify rj' physics in ways relevant for this experiment [T^ and boost 
1]' Dalitz to a more important level. Direct decay of vector mesons p° and u will also be 
a strong source of pairs but fairly near their free space masses. Production of a nonzero 
four volume of hot matter generates contribution from pion annihilation tt^tt^ e~^e~, 
whose strength depends on the average temperature. In this four volume there are myriad 
hadronic reactions able to produce pairs. A systematic treatment of decays and annihilation, 
i.e. lowest order treatment in the strong coupling has been done |T^. Going beyond lowest 
order is required for this discussion as np ire^e^ comes from a two-loop contribution. 
Elastic scattering of vr's with p's is dominated by resonant ai formation as the average ^/s 



for T ~ 150 MeV is near the centroid of the rather broad oi mass distribution. Since 
the ai has a radiative decay channel as well, it follows that scattering can produce real or 
virtual photons. For real energetic photon production, the reaction np —>■ ai ^717 has been 
shown to be very important |T^JT5[] . As ^ phase space shrinks to a point and this 
mechanism contributes nothing. Zero invariant mass for virtual photons does not require 
zero energy as and g ^ could separately be very large. Consequently the amplitude for 
TTp Txe^e'^ will go like for small mass. Near the rho, vector dominance takes over 

and the amplitude goes like 1/(M^ — + inipT p). A completely equivalent and alternative 
approach to estimating the role of the ai as considered here would be to compute ai — > 7re^e~ 
starting from a thermal distribution of charged Oi's properly modified by a Breit-Wigner 
distribution in mass [1^ owing to its finite width. 



More general techniques of field theory at finite temperature unambiguously determine 
the total contribution The rate for producing electron pairs is related to the imagi- 

nary part of the retarted photon self-energy which, to one-loop order, corresponds to decay 
and annihilation. The pion-resonance scattering process being of two-loop order would re- 
quire additional renormalization and the added nuisance of an additional Matsubara sum, 
but the prescription is clear. As the field theoretic methods are completely equivalent to 
kinetic theory at given order in the coupling, the latter is chosen for convenience. It builds 
upon evaluating tree-level Feynman diagrams and folding thermal momentum distributions 
with transition amplitudes to arrive at an average rate in medium for a given process. For 
instance, the dominant ai resonant graph contributing to 7rp vre+e^ is shown in Fig. [l|. 
Interference and other meson-exchange effects will be reported upon separately ||2^. The 



appearance of the intermediate rho suggests using a vector-dominance form factor whose 
width is taken to be a function of temperature to include collision broadening. It is com- 
puted just as in Ref. |]19| but including finite pion chemical potential and then parametrized 



for 100 < T < 200 MeV by 

Tp{T) = Tp+\a + hT + cT^] (1) 
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where Tp is the free space width and a = 0.50 GeV, b = —7.16 and c = 30.16 GeV ^. 

An effective lagrangian for axial-vector-vector-pseudoscalar interaction is chosen for sim- 



phcity to be the following [jn| 



LavcP = 9Av<t>Af, [{p^ ■ Pv)g^'' - p'^^Pv] K (2) 

where A, V, and are respectively, axial-vector, vector and pseudoscalar fields. A treatment 
of this interaction from a chiral lagrangian approach has yielded results for photon produc- 
tion roughly consistent with those from this less complicated parametrization |I^p!5|. The 



coupling constant is adjusted to give a decay Ta^^-Kp = 400 MeV pT| to match the Particle 



Data Group value and the coupling of ai to vr and 7 is then taken to be the vector- 
dominance value of gaiTTpi^/ fp), where fp is the coupling of rho to pions. Numerically these 
turn out to be gamp = 16.1 GeV"^ and gaiw'y = 0.81 GeV~^ for rria^ = 1230 MeV. The result- 
ing prediction for the radiative decay width Ta-^^^T^y is 1.9 MeV, which is somewhat larger 
than the 0.640± 0.246 MeV obtained from a sole measurement p3|. Consequently, rates for 
photon or dilepton production from tt + p scattering through the ai might be corresponingly 
lower than this approach would indicate. 

The invariant rate for thermally producing a lepton pair of mass M and individual three 
momenta and via the process p{pa) + T^iPb) —* '^{pi) + ^+ ^- can be written generally 
as 

(Pa +Pb-Pi-U- f--) S + i-Y - M^] (3) 

where AT is an overall degeneracy factor, dua = d^pa/ [{271^2 Ea] and so on, / is the Bose- 
Einstein distribution, / = 1 + / to account for medium (Bose) enhancements and | 7W p is 
the initial spin averaged and final spin summed squared matrix element. Integration over 
the full individual lepton three momenta can be additionally performed to arrive at the total 
rate for given mass. However, the differential rate is more useful here since limited trans- 
verse momentum and rapidity ranges can then be integrated to approximate experimental 
configurations. 



If we ignore the Bose enhancement then the rate can be simphfied as 

da 



(4) 



dM^ dH+ dH. 

with z = \fsjT and K\ is the first modified Bessel function. For completeness, the cross 



section in brackets is 



^^^-dM^H^dH. ~ 8(2vr)5A' /2(J,m2,mg)^ (^^» + - ^+ - " "^?) 



x5 



{l+ + i^f-M^ . (5) 



First, we integrate Eq. @ to get the total rate dN / d'^xdM'^ . To faciliate comparison with 
other rates from the literature we temporarily ignore collision broadening and set the pion 
chemical potential to zero. In Fig. |^ the rates from txtx annihilation are compared with 
Tip — > 7re^e~. Formulas for annihilation are not shown because they are by now textbook 



expressions and have been recently used in a similar study For additional com- 
parison we show also an estimate of pion scattering with bremsstrahlung All three 
contributions become roughly the same around 350 MeV mass, which is the lower end of 
the window of the observed excess over hadronic decays. The mechanisms have different 
Pt distributions for the leptons which will become important when kinematic restrictions 
are imposed. At the rho mass the two-pion annihilation is a factor of four above the vrp 
scattering. We can expect this number since it is precisely the ratio of thermal production 







rates of neutral rho mesons via the two mechanisms tt+tt p° and ai ixp 

dN^^^P° IdN'^^^^p'' 

at T = 160 MeV. 

In order to compare with experiment an integration of these rates over a spacetime 
evolution must be preformed. We assume for utmost simplicity boost-invariant expansion 
p6| and neglect any transverse flow effects. This will provide a first estimate only while the 
need for dynamic model calculations is duely noted. Choice of specific times or temperatures 
must be made. Here we take % = 188 MeV, = 160 MeV, Tf = 140 MeV and /i^ = 135 

6 



MeV. The rho is again assigned a collision broadened width. Total yield is the sum of 
contributions from the mixed plus cooling phases p|. Written out in detail, it is 



dydM 2 \TJ ' ' ' J E+ E 



77. 77. (T = T ] 

+ 'd^xdMdH+dH-^ 



dT d^^ dH^ 



E, E^ 



d^xdMdH^dH- 



(7) 



where r is the ratio of degrees of freedom in QGP phase to that of hadron phase and is 
~ 12. The pseudorapidity density will be assumed to be equal to the rapidity density 
dN/dydM ^ dN/drjdM. 

Attempting to resemble the experiment as closely as possible lepton transverse momen- 
tum integrations are limited to pt > 200 MeV/c and pseuodrapidities by 2.1 < t] < 2.65. 
There is also a cut in the experiment on angular separation for the pairs O > 35 mrad 
to avoid resolution difficulties. We ignore this limitation here as it should only have small 
influence on the result for masses greater than 200 MeV. Finally, we normalize just as in the 
experiment by dividing the yield d'^N/drjdM by the average charged particle pseudorapidity 
density (rfA^charged/'^'?) = 125. 

Having established the ranges kinematics "visible" by the experiment, the rate from 
Eq. can be obtained and then the convolution of Eq. performed. The resulting nor- 
malized yield after performing a mass-resolution smearing might then be directly compared 
with the distribution observed by the CERES collaboration. Partial and total yield esti- 
mates are shown in Fig. |^ as dashed (vrvr annihilation), dot-dashed [irp scattering), dotted 
(hadronic decays as presented in Ref. 0) and solid (total of annihilation, scattering and 
decays). Inclusion of the Oi seems to account for part of the excess but still leaves the door 
open for other possibilities. Additional pairs could be coming from tttt — > p7* and tttt iT'-f* 
p7| . They could also be coming from decays not already considered and possibly even a 
bremsstrahlung component pB|. 



As a brief summary, we have pointed out that the recently observed dielectron signal 
being in excess of hadronic decays in S+Au collisions at 200 GeV/u is likely coming not 



only from tttt annihilation as was suggested, but could also be partly due to pion plus 
resonance vrp Tre'^e" scattering. The picture emerging for thermal production of low- 
mass dielectrons becomes richer and more challenging as it calls for next-to-leading order 
treatments of the photon self-energy, i.e. 2^3 body reactions (including the leptons). This 
is not too surprising however, as there are many hadronic reactions of type a + b ^ c + i'^i" 
for the full set of light mesons. The present result for irp scattering suggests that moving 
towards quantitative interpretation of the CERES data will require consideration of such 
processes. 
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FIGURES 

FIG. 1. Dominant resonant contribution to irp ire^e" scattering. 

FIG. 2. Thermal rate for producing electron pairs at T = 160 MeV for the process 
(ttp — iTe'^e~) in the solid curve, pion bremsstrahlung in dotted curve and tttt annihilation in 
the dashed curve. 

FIG. 3. Normalized lepton pair yields observed in experiment by the CERES collaboration 
as compared with pion+resonance scattering (tt + p — > 7re"'"e~) in the dot-dashed distribution, 
estimated hadronic decays through event generation taken from Ref. [6] (dotted), tttt annihilation 
in the (dashed) and with the sum of annihilation, scattering and decays (solid). 
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Figure 1 



K. Haglin, ' 'Excess electron pairs from heavy-ion collisions at CERN 
and a more complete picture of thermal production." 



K. Haglin, ' 'Excess electron paks from heavy-ion collisions at CERN 
and a more complete picture of thermal production." 
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K. Haglin, ' 'Excess electron paks from heavy-ion collisions at CERN 
and a more complete picture of thermal production." 
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